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Abstract 

In this paper we check the validity of the ideas of double scaling as given by Ball and Forte in a 
semi inclusive process like J/tp production at HERA, in different kinematical regions, for low values of 
the Bjorken variable x. In particular, we study J/ip production in the inelastic and diffractive (elastic) 
regimes using the double scaling form of the gluon distribution functions. We compare these predictions 
with data (wherever available) and with other standard parameterisations. We find that double scaling 
holds in the inelastic regime over a larger kinematic region than that given by the analysis of the proton 
structure function FJf. However, in the diffractive region, double scaling seems to suggest an admixture 
of hard pomeron boundary conditions for the gluon distribution, while predicting a steeper rise in the 
cross section than suggested by present data. 
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1 Introduction 


The systematic study of the dynamics of perturbative QCD at low values of the Bjorken variable x has 
thrown up many new ways of analysing the structure of the proton through a measurement of the structure 
function F %. In particular, one of the methods that provides a clean study of the scaling properties of 
structure functions as Q 2 (the virtuality of the photon) increases and x decreases is the so-called double 
scaling analysis of Ball and Forte JjJ. In their work, in the kinematic region of very large values of the 
photon-proton centre of mass energy squared, S , the Altarelli-Parisi (DGLAP) evolution equation were 
solved and it was shown that the gluon distribution (which is the dominant partonic distribution in the 
region of interest) shows double scaling in terms of the variables 


„ io, t /in xq/x 


a) 


where the starting scale for Qq in to = In -0- can be just a little more than Q q = lGeV 2 , A = A qcd, 
Xq = 0.1, typically. This asymptotic analysis is a straightforward consequence of perturbative QCD 
evolution of structure functions at low x , and provides excellent agreement with the data for Q 2 > 5 GeV 2 . 
However, taking into account the fact that other analyses, the most common being the BFKL approach Q, 
also provide a mechanism for predicting the behaviour of F 2 at low x , a case has been made for studying 
the validity of some of these different approaches in processes that are not as inclusive as F 2 . One such 
process where the validity of the double scaling hypothesis has not been tested is J/production. In this 
paper, we propose to remedy this by studying J/\P production at high centre of mass energies of the kind 
available at HERA, using the double scaling form for the gluon distribution function. 

The issue of J/'F production at large centre of mass energies available at HERA (y/S ~ 200GeV) has 
become an important aspect of QCD because, generally, production of heavy quarkonium states in high- 
energy collisions is an important tool for the study of its various perturbative and non-perturbative aspects. 
The production of heavy quarks in a hard scattering process can be calculated perturbatively. However 
the subsequent creation of a physical bound state (J/\F or other quarkonia) involves non-perturbative 
QCD. This latter process is treated in a factorised approach where the short and long distance physics are 
separated out. In particular, the evolution of a quark-antiquark pair into a physical quarkonium state is 
described by the color singlet or color octet models to which we will return in more detail later. 

In this paper, we shall address the issue of J/\F production at low values of the Bjorken variable x. At 
low values of x, the gluon distribution plays a dominant role in the determination of the cross section as 
mentioned earlier. In particular, we concentrate on the specific kinematic region of very large values of 
the photon-proton center of mass energies, y/S, in e-p collisions. We calculate the J/’F production cross 
sections at low x, using the asymptotic form for the gluon distribution function, given by the (leading order) 
double scaling analysis in different kinematical regimes, namely for inelastic and diffractive scattering. Till 
now, most calculations for J/’F production use the standard parametrisations for the gluon distributions, 
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based on a fit to the known data. The advantage of using the double scaling form is that it is, in a certain 
sense, derived from QCD and has no fitting parameters other than the starting values of Q 2 and x from 
where the distributions are evolved. We have, in fact, checked the dependence of our predictions on these 
starting values in order to be able to quantify the degree of uncertainty that arises from these parameters, 
for diffractive J/T photoproduction. 

The paper is organized as follows. In section 2, we briefly summarize the phenomenon of double 
asymptotic scaling of the gluon distribution function. In section 3, we present our results for J/T production 
in the inelastic regime. We then, in the next section move on to the elastic regime, where BFKL dynamics 
is supposed to be dominant and make our predictions regarding the cross sections. Although it may 
seem inconsistent to use the gluon distribution function obtained by solving the Altarelli-Parisi evolution 
equation in the double scaling limit in the diffractive regime where BFKL dynamics takes over, it may be 
pointed out that these distribution functions are merely used as parametrizations of the data. 

In the last section we discuss our results and various open problems. 

All the experimental data included in this paper have been obtained from |j) except where otherwise 
stated. 


2 Gluon Distribution at low x 


There have been two distinct approaches to the study of the HERA data, the ’standard’ Altarelli Parisi 
(DGLAP) [|J evolution equation approach and the attempt to reconcile the data with BFKL dynamics 
||]. In order to study the effects of summing a s In A unaccompanied by InQ 2 , the preferred approach has 
been to study the BFKL equation which generates a singular x~ x behaviour for the unintegrated gluon 
distribution f(x , kj.), with A = a s 41n2 (for fixed, not running a s ). For running a s , A ~ 0.5. Q. From this 
an asymptotic form for Fi(i = 2, L) can be inferred. We shall come back to the issue of BFKL dynamics in 
more detail when we deal with diffractive J/T photoproduction, later in the paper. In the other approach 
one attempts to describe the data through the DGLAP evolution equation to the next-to-leading order 
approximation. Here again, the data imply a steep gluon distribution with the gluon density rising sharply 
as x decreases, even for comparatively low values of Q 2 . Ball and Forte Q have used this approach to 
exhibit the scaling properties of the gluon distribution function and hence the proton structure function 
F'i at low x, generated by QCD effects and have shown that the HERA measurement of F 2 (x, Q 2 ) is well 
explained by this approach. In Jlj] it was shown that the solution of the traditional DGLAP equations 
in the low x regime exhibits double scaling in terms of the variables defined in (1). Double asymptotic 
scaling results from the use of the operator product expansion and the renormalisation group at leading 
(and next-to-leading) order, and predicts the rise of F 2 on the basis of purely perturbative QCD evolution. 
The asymptotic behaviour of the gluon density in this approach is given (to leading order) by 
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which leads to the asymptotic behaviour of F 2 : 
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n c and n f being the number of colors and the number of flavors respectively. N is a constant that depends 
on the input gluon distribution. In our analysis, we have taken N = 3.24, corresponding to soft boundary 
conditions on the gluon jlj as favoured by HERA data j5| , in the region of x that we are interested in, for 
inelastic J/\k production. For J/\R photoproduction in the diffractive scattering regime, we shall see that 
harder boundary conditions are preferred and in this kinematical region, we have taken an intermediate 
pomeron, corresponding to N = 1.45. 

In our calculation, for the J/\R production cross sections in this double scaling limit, at small values 
of x, we have used the values Q 2 = 1.12GeF 2 ,A = 0.248GeF and the starting value of Xq = 0.1. For 
the values of the variables x and Q 2 for which a 1 > 1 and p 2 > 2, double scaling of the proton structure 
function F 2 is confirmed to a fair degree. However, in the semi-inclusive J/'k production process, we shall 
not restrict ourselves to these kinematical regimes, but rather try to check the region of validity for double 
scaling where data is available. 


3 The Inelastic Region 


In the color singlet model approach, the dominant contribution to the quarkonium production comes from 
quark-antiquark or gluon-gluon fusion. The fusion process produces a heavy quark pair which forms a 
physical color singlet state with an appropriate color singlet projection. This fusion contribution to J/*k 
production in the color singlet model has been computed in Q. It takes place via the subprocess 7 g —► J/'S’g 
where the gluon in the final state of the partonic scattering process carries off the color charge. The final 
cc state is a color singlet and the cross section for the process 7 N —> J/'SfX can be calculated from the 
above subprocess and is given by 
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Figure 1: Color Singlet contribution to inelastic J/\& production. 


with s + t + u = m 2 j , and Y J ee the electronic width of the J/’F which is related to its orbital wave 
function at the origin (Fg e = 5.26KeV in our calculations). Here, s, t , u are the usual Mandelstam variables 
for the partonic subprocess and s is related to the total centre of mass energy S by s = xS ,x being 
the fraction of the incident nucleon momentum carried by the gluon. It is more convenient to re-express 
this differential cross section in terms of the variables 2 the elasticity parameter, and P 2 , the transverse 
J/\F momentum squared. The elasticity parameter 2 is defined by 

P JN ,.Pn 
Z ~ Q.P n 


where Q = P 7 = momentum of the incident virtual photon. 2 denotes the fraction of the energy 
transferred to the J/T in the process and in particular in fixed target experiments, reduces to 2 = Ej/^/E^ 
In terms of 2 and Pt , we may write the Bjorken variable as 
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x = 


m " 
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It is more transparent to reexpress this differential cross section in terms of the variables 2 and P^. In 
terms of these variables, the differential cross section is given by, 

d 2 <j xg(x)z(l — z)nijA - 


dP 2 dz 


[mj(l — 2 ) + Py] 2 
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This equation can be used to evaluate the differential cross sections jpr, jp and also the Pp integrated 
total cross section. It may be remarked here that the above expressions show that the cross sections are 
finite for non zero Pp for all values of z, and for Pp = 0, it diverges in the limit when z = 1. However, 
this value of z is excluded from the domain of our interest, as there is a strict upper bound on z for this 
inelastic process, viz. the restriction that the momentum transferred in the partonic subprocess be above 
the QCD scale so that perturbative QCD is valid; this restriction, \t\ > Q § gives 


(1 - z) > 


Ql-P$ 


( 8 ) 


Ql+m 2 j 

and hence, corresponding to Qq ~ 1.12GeV 2 , we have z < 0.90. 

Before we present the results of our calculation, a word about the applicable region of phase space for 
the double scaling form of the gluon distribution function is in order. We have already remarked that 
the predictions of the proton structure function F% at HERA shows that asymptotic behaviour sets in, in 
terms of the variables er and p defined in the beginning, for er 2 > 1 and p 2 > 2. As can be seen from 
the expression for x, imposing the cuts z < 0.9, we cannot really go to sufficiently small values of x for 
large values of Pp, which is essential for the validity of the asymptotic form of the gluon distribution. To 
be well within the double scaling region in x and Q 2 , for such large values of Pp, one thus has to impose 
much stronger cuts on z. However, we shall not use such cuts in our analysis as we try to go beyond the 
limits for a and p as obtained from the analysis of F%. 

We have analysed the behaviour of the total cross section and the Pt and z differential cross sections 
for J/T production. The results of the analysis are shown in Figs.(2),(3),(4). In our calculation, we have 
taken m c = 1.4 and the starting scale Xq is taken to be 0.1. For the calculation of the total cross section, 
we have imposed the cuts 0.1 < z < 0.8 and Pp > lGeV 2 , whereas for the Pt differential cross section, the 
cut on the elasticity paramter has been imposed as 0.1 < z < 0.9 in order to match with experimental data 
from Q. Cuts on Pt are generally imposed as in the region of small Pt, the non leading order calculations 
are not fully under control. Also, the lower bound on z is required in order to suppress contribution to 
,1 /\H production arising out of B decays. We have checked that on inclusion of a K factor of 1.7 in the total 
cross section data (as in ||), the results are in fair agreement with extrapolated data points ||. From the 
results it is also clear the one can extend the double scaling validity region to kinematic ranges that go 
beyond that obtained from the analysis of the proton structure function F%. Indeed, one gets fairly good 
agreement with the data for values of the variable a 2 > 0.4, p 2 being typically greater than 2. 

We should mention here that in the expression for the total cross section and dcr/dz, one has to in 
principle integrate over all values of Pp > 1 GeV 2 . From the expression for x, it is clear that to remain at 
low values of x , one cannot really go to large values of Pt {Pt > 6 GeV, say), without severely constraining 
the inelasticity parameter z, as we had pointed out. However, as can be seen from the expression for the 
cross section, in the approximation that Pp is large compared to m 2 , the denominator falls off as Pp, and 
hence the contributions from the large Pp region is negligible for such large values of Pp. Hence, in our 
numerical calculations, we have restricted ourselves to Pt < QGeV. 
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Figure 2: Total cross section for inelastic .I/T production. The solid line corresponds to the Ball and 
Forte gluon distribution, the dashed line to the GRV and the dotted line to the CTEQ parametrisations. 
The dashed-dotted line shows a fit to the data with the Ball and Forte distribution assuming a K-factor of 
1.7. 


The z and distributions are shown in Figs(3) and (4). From the P%. distribution, it is seen that the 
double scaled form of the gluon distribution function underestimates the data to some extent, however, 
the K factor (enhancement factor for the cross sections on inclusion of the NLO corrections) might take 
this into account at least in part. In the calculation of the differential cross sections also, we have, 
as mentioned earlier restricted ourselves to the region Pf > lGeV 2 as in the limit of extremely small 
transverse momentum, there are large negative contributions coming from NLO corrections and the cross 
section falls sharply as Pr —> 0. 

In our calculations, we have not considered ” color octet” contributions to the J/T cross sections. In the 
color singlet model, that we are considering, one conventionally assumes that the color singlet cc state is 
produced in the 3 Si state initially and then evolves into the physical J/ , F state. This model however grossly 
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Figure 3: z Differential cross section for inelastic J/*F production at \JS lp = lOOGeE The solid line 
corresponds to the asymptotic gluon distribution, the dashed line to the GRV fit and the dotted line to 
the CTEQ parametrisation. 


underestimates the J/*F production data at the Tevatron jfj] . To reconcile with the data, a new mechanism 
for J/\F production was suggested in 0 , the color octet mechanism. Whereas in essence the ’’color 
singlet” model ignores the relative velocity v between the quarks in the bound state, this quantity is not 
often negligible and may give rise to considerable amount of corrections to the results. In [|ll| an expansion 
of the quarkonium wave function in powers of v and corrections involving the octet components of the wave 
function have been carried out. However, as has been pointed out in []l2| color octet contributions may 
become important only in the diffractive limit, i.e z close to 1 , for small Pt- It was further demonstrated 
in ]l2j that at HERA experiments, in the inelastic region, i.e for z < 0.9, P^ > 1, the data can be well 
accounted for by inclusion of the next to leading order corrections to the cross sections, and hence the color 
octet contributions in these energy ranges do not seem to play a very significant role. Comparison of the 
color singlet and color octet contributions to r j^ with experimental data has been carried out in (l^] and 
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Figure 4: P^ Differential cross section for inelastic J/\F production at \/S lp = lOOGeR. The solid line 
corresponds to the Ball Forte distribution, the dotted line to the GRV and the dashed line to the CTEQ 
parametrisations. 


this analysis shows that the J/\F spectrum is well accounted for by the color singlet model alone and the 
octet model shows a marked increase in cross section at large values of z, which is clearly in conflict with 
experimental data. Hence we are justified in dealing only with the color singlet model in our calculations. 

The choice of scale as Q = Pt/2 would have in part accounted for the K factor mentioned earlier; 
however, with this choice of scale, small values of Pt in the double scaling regime cannot be reached, as 
from the form of the gluon distribution function (2) defined in terms of the variables er and p in (1), it 
is clear that it diverges for Q = Qo- We have thus chosen our scale as the factorization scale, namely 
Q 2 = 2 m 2 c . 

Also, at HERA energies, contributions to the J/\F production cross section might arise from B decays, 
however these start contributing significantly at extremely low values of z, namely z < 0.1 ||l^| We have 
ignored such small values of 2 in our analysis, hence this effect will have negligible effect on our results. 








Before ending this section, we would like to point out that although we have considered here only the 
fusion process contribution to J/\I/ production, at large values of Pt, the fragmentation process, i.e the 
fragmentation of gluons and charm quarks, also become an important source of J/'f' production. Even 
though these contributions appear at higher orders of the strong coupling as compared to the fusion con¬ 
tributions, they become important at large values of the J/H/ transverse momenta due to an enhancement 
factor of Pf, /rrr c . m c being the charm quark mass. Consequently, fragmentation process contributions 
may dominate over fusion at comparatively large values of Pt- Fragmentation functions are computed 
perturbatively, at the initial scale of the order of the charm quark mass and the Altarelli-Parisi equations 
are then used to resum large logarithms in Pt/to c . Gluon and charm fragmentation functions have been 
calculated |Q,||l^] and it has been shown [7i] that at HERA experiments, the charm quark fragmentation 
process is expected to dominate over the quark gluon fusion process at sufficiently values of Pt > 7 GeV. 
Calculation of the fragmentation functions in the double scaling regime remains an open issue and we hope 
to address this elsewhere. 


4 The Elastic Region 



Having looked at J/\I/ production in the inelastic limit, we now move over to the diffractive region, i.e 
for values of the elasticity parameter 2 close to 1. The description of the scattering process via the standard 
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photon gluon fusion mechanism breaks down in this limit, as the elasticity parameter takes values close to 
1, and the formulae of the previous section are no longer applicable. In this kinematic regime of diffractive 
scattering, BFKL dynamics takes over, and in particular, one has to sum over the gluon ladder diagrams 
as shown in Fig(5). In this elastic region, the scattering amplitude occurs in a factorized form, 
which is essentially due to the fact that in the proton rest frame, the formation time of the J/\F is much 
greater than the time of interaction of the photon with the proton. The photoproduction amplitude with 
the incident photon almost on shell is described by a two gluon exchange scattering process where the 
transverse momentum transfer t ss 0 for the diffractive process. The amplitude for the process in [Fig. 5] 
was evaluated in the leading log approximation in perturbative QCD in jl(j] and depends on the square 
of the gluon distribution function. Subsequently, this formula has been improved upon in ||l7f where 
the transverse momentum of the gluon ladder is included in the analysis, thus incorporating BFKL [^j] 
dynamics. Dependence on the square of the gluon distribution function makes this cross section extremely 
sensitive on the latter and a precise measurement of the cross section will shed light on the behaviour of 
the gluon at very small values of x. It may be mentioned that here we are using the DGLAP evolved gluons 
in the double scaling limit as a mere parametrization of the data, and hence it is not inappropriate to use 
them in a regime where BFKL dynamics plays an important role. In ]l6| the amplitude for the Feynman 
graphs of [Fig. 5] were evaluated and found to be 


dcr(7P —> J/’F + p) 
dt 


{t) 192a 


- ^2Q 2 - |Qt| 2 

xg{x,Q ) 


{Q 


2\2 


(9) 


where a, a s and denote the electromagnetic coupling, the strong coupling and the electronic width 
of the J/’F as before and the quantity F denotes the two gluon form factor that is apriori unknown. 
To the first approximation it can be treated as the electromagnetic form factor, but its nature is strictly 
determined experimentally. We shall however restrict ourselves to the domain of photoproduction of J/\F in 
the diffractive regime, i.e where the photon is almost on shell and the momentum transferred, t is extremely 
small, t ~ 0. At such vanishingly small values of t, this form factor can be approximated to unity |lfj. 
Here, Qt denotes the transverse momentum of the J/\F particle, which is typically small, ~ 0.3GeV\ The 
quantities x and Q 2 is defined as 


Q 2 


- (|Q“| + TOj) 



( 10 ) 


where S is the photon-proton centre of mass energy. As we wish to consider the photoproduction limit, 
we set the photon virtuality to be extremely small, and in our calculations, we have chosen the scale 
Q 2 = 2AGeV 2 In this expression, only longitudinally polarised gluons have been taken into account as 
these are the only ones that contribute in the extremely high c.nr energy limit. We note, first of all, that 
the dependence of the cross section on the barred quantities ensure that perturbative QCD can be applied 
even in the photoproduction limit, due to the heavy mass of the J/*F . Secondly, the expression for x tells 
us that at sufficiently large values of S, the extremely small limit of x can be probed. 
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This formula (9), however has been derived in the leading In Q 2 approximation, without the inclusion 
of the transverse momenta of the exchanged gluons. The analysis of the above process including the gluon 
transverse momenta has been systematically carried out in [jlT]]. In their calculation, the cross section 
is expressed in terms of the unintegrated gluon distribution, f bfkl{x, fc^), which satisfies the BFKL 
equation, and is related to the conventional integrated gluon density by the relation 

rQ 2 


r v jia 

xg(x,Q 2 ) = J -gffBFKL{x,kT) 


Inclusion of the gluon fcf, modifies the scattering cross section to 


dcr( 7 P —> J/\F + p) 
dt 


167T 


with, 


and B 2 = 


2 „ f dkj, 

A = i2ir mja s B J 


k t 


(Q 2 )(Q 2 + k 2 ) 


dxg(x , kjd) 


dk t 


( 11 ) 


( 12 ) 


(13) 


The amplitude A can be written in terms of the conventional gluon distribution as 

,, o 2 D l ~ x 9( x ’Ql) , f°° dk T dxg(x,k%)~\ 

A = .2, m ja ,B (- 53 — + j Qi ((5W + ^ ) d e T j (n> 

The lower cutoff Ql has to be put on the integration due to the undetermined form of the gluon distribution 
function at vanishingly small values of fcf.. This form of A can be plugged in (12) in order to obtain the 
total cross section. 

We have obtained the values for the total cross section for diffractive photoproduction of J/T mesons, 
using (12) and (14). From our analysis, it seems that using the soft pomeron normalisation in this region 
a la Ball and Forte jl|, the double scaling analysis grossly overestimates presently available data, by a 
factor of about 5. Hence, we conclude that in this region of extremely smally values of x, x < 10~ 3 , there 
is some admixture of the hard pomeron boundary condition, and consequently to analyse this region, we 
have chosen a normalisation corresponding to some intermediate pomeron, namely, N = 1.45. 

As a lower cutoff for the transverse momenta of the gluon distribution function, we have used the value 
Q 2 l = 1.5GeH 2 . This value of Q 2 L has been chosen so that it is just above the starting scale for the gluon 
distribution (namely Q g = 1.12GeH 2 ) in the double scaling analysis. We will come back to the issue of 
variation of the cross section with this cutoff in a moment. Note that in our calculation, we have considered 
the forward scattering amplitude, namely t = 0. There may be corrections coming from consideration of 
a non zero value of the momentum transfer, however, these corrections are expected to be small in 


the region of x that we are interested in, namely x < 0.01 117 . One might also consider the corrections 
arising out of the relativistic effects of the J/'F wave function, those coming from treating the quarks in the 
J/’F bound state non relativistically, i.e assuming mj = 2m c , where m c is the mass of the charm quark, 
and also those arising out of radiative corrections to the lowest order process considered here. Analysis 


of these correction factors have been carried out in 117 and it has been shown that whereas the radiative 
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Figure 6: Total cross section for diffractive J/\F photoproduction. The lower, middle and upper lines 
corresponds to x$ = 0.1, 0.15 and 0.2 respectively. 


correction is negligible for soft gluon emission, as is the case here, the total correction factor arising from 
relativistic effects in the J/*F wave function and the motion of the quarks inside the charmonium bound 
state is of the order of unity. 

The results of our analysis is shown in Fig [6]. From the available data points ||, we find that the 
double scaling form of the gluon distribution function predicts a slightly steeper rise in the elastic cross 
section than favoured by the experimental data. However, experimental values for larger centre of mass 
energies than are currently probed may throw light on the scenario for diffractive J/\F photoproduction. 
We have also plotted the variation of the cross section with the choice of the starting scale xo in the gluon 
distribution function in the double scaling limit. It can be seen that change of Xq causes a considerable 
amount of fluctuation in the J/\F production cross-section. We have checked that there is similar fluctuation 
in the cross section as one varies the starting scale for the gluon distribution function, namely Qg, though 
the results are not explicitly indicated here. We have also explored the sensitivity of the cross section to 
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the lower cutoff in the transverse momentum (Ql) used to compute it. It is found that the cross section 
is highly sensitive to fluctuations in Q 2 L , varying from about 70% at V~S = lOOGeK to nearly 75% at 
\fS = 200GeK. This is in contrast to the smaller variation in cross section noted in jf(| using the GRV 
and MRS gluons, as Q 2 L is varied. 

5 Conclusions 

In this paper, we have studied the validity of the double scaling ideas to a specific semi-inclusive process 
like J/T production. Our results show that for the inelastic region, the data is fairly well accounted for 
by the double scaling form of the gluon distribution function, even in a region beyond that confirmed by 
measurement of the proton structure function F 2 . However, in the diffractive scattering regime, double 
asymptotic scaling while predicting a steeper increase in the cross section than from present data indicates 
the admixture of hard pomeron boundary conditions (as opposed to the soft boundary conditions |(J 
that we assumed while analysing the inelastic scattering data). As such, it is not very clear, at least 
from the present analysis as to which behaviour takes over at extremely small values of x. and diffractive 
J/\k photoproduction may provide some crucial hints, as it is extremely sensitive to the gluon distribution 
function. 

To our knowledge, this is the first attempt to apply the idea of double scaling to a semi-inclusive 
process. The advantage of using double scaling, as pointed out in the introduction is that it does not 
need any given form form for the starting gluon distribution g(xo,Qo) as with other parametrisations. 
We have studied whether the region of validity of double scaling in a totally inclusive process like the F 2 
measurement changes when a semi-inclusive process like J/'F photoproduction is analysed. In the case of 
diffractive scattering, we have also used the data to constrain the starting value of x, i.e Xq and found 
significant change in the value of the total cross section, a similar change being present if one varies the 
starting scale for the gluon distribution function, i.e Qq. 

There have been many suggestions, including a very elegant one by Mueller jl^] for studying processes 
that would clearly distinguish QCD evolution through DGLAP and some other mechanism like BFKL. 
The basic thrust of all these suggestions has been to study semi-inclusive processes and we have taken a 
first step in this direction by studying the validity of double scaling in J/T photoproduction where data 
with reasonably good statistics is already available from HERA. It would be useful to make a systematic 
study of all such semi-inclusive processes for which data exists, in different resummation schemes (DGLAP, 
BFKL, CCFM, etc) in order to look for the validity of these different methods. This would then give us 
a qualitative understanding of the physics of the different processes and consequently that of the various 
resummation schemes. These are some of the issues for future study. 
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